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We use the Density Matrix Renormalization Group technique to compute cor- 
relation functions of the Za-chiral Potts quantum chain in the massive regimes. 
Chains of up to 70 sites are used and a clear oscillatory behavior is found. We 
also check the relation between the density matrix spectrum and the Corner 
Transfer Matrix spectrum pointed out by Nishino and Okunishi, using Baxter's 

ly^ ■ low-temperature formulae for the chiral Potts Corner Transfer Matrix eigenval- 

lO ' ues. We find very good agreement. 

(N' 

^ ! 1 Introduction 

^ ■ We report numerical calculations which apply the density-matrix renormal- 

c^ ■ ization group method (DMRG) for the first time to the chiral Potts quantum 

chain. Although the chiral Potts model is integrable, various interesting prop- 

_^ , erties are not yet accessible to analytical calculation and one has to resort to 

Ch ! numerical calculations in order to get more complete information. We shall 

^ ■ concentrate on two aspects: first, we show clearly that in the high-temperature 

massive phase the correlation functions are oscillating. Second, we check the 

. relation p ~ ABCD between the density matrix and the Corner Transfer 

rN ; Matrices (CTM), pointed out by Nishino and Okunishi |30|. For the CTM 



_c^ \ eigenvalues of the chiral Potts model we use Baxter's low-temperature expan- 



sion formulae ||T0|| . 



2 The integrable chiral Potts model 

The Zjv-symmetrical chiral Potts model has been introduced by Ostlund |^ 



in order to describe incommensurate phase transitions in surface layers, e.g. 



Krypton monolayers on graphite. Centen et al. ||T5[ introduced and analysed 
the chiral Potts quantum chain. Very soon it was realized that this class of 
models can be reformulated such that very interesting mathematical structures 



emerge: Howes, Kadanoff and den Nijs |^ modified the A^ = 3 quantum chain 
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to make it self-dual and discovered that in their version for a particular chi- 
ral angle the lowest gap depends linearly on the inverse temperature. Then 
von Gehlen and Rittenberg |2^ constructed a whole set of Zat symmetrical 
quantum chains which for any A^ satisfy the Dolan-Grady [|1^ integrability 
conditions. These quantum chains are now usually called "superintegrable". 



Perk noticed [33] that the Dolan-Grady conditions are just conditions for gen- 
erating the algebra which Onsager in 1944 had introduced to obtain the free 
energy of the Ising model ||31 . 



Onsager's 1944 technique for solving the Ising model made use of specific 
Z2-features and it has not been possible to adapt it for solving the Z^v- 
model. The actual solution for the eigenvalues of the superintegrable model 
was achieved differently: In 1987 Au-Yang et al. invented a two-dimensional 
Zjv-symmetrical Yang-Baxter-integrable statistical model, which contains the 
chiral Potts chain hamiltonian as the derivative of its transfer matrix. Inver- 
sion relations and functional equations were found for this 2-dim. model and 
were used to obtain the ground state energy and also analytic formulae 
for the excited levels |1|, |T6|. Hints that the chiral Potts model is related to 
the six-vertex model at quantum group parameter at root of unity came from 
Korepanov |^ and the connection was shown by Bazhanov and Stroganov[^ 



We first summarize some basic notions of the 2-dimensional integrable chi- 
ral Potts model of |]^. Fig.l shows part of a diagonal square lattice and rapidity 
lines used to define the Boltzmann weights Wp^q{aj — cr'j) and Wp^g{aj — cr'j+i) 



q - 




Figure 1: Diagonally drawn square lattice with Zjv-spins aj on the lattice sites. 
Dashed lines: p- and g-rapidity lines. Boltzmann weights Wp^g (^aj — a'j) act 
in SW-NE-direction, Wp,, {a^ - a'j+i) in SE-NW-direction. 

which depend on both rapidities via six functions Xp, yp, fip, Xg, yg, fig (for brevi- 
ty, we write x{q) as Xg, etc.): Wp^g{0) = Wp^g{0) = 1 and for n = 1, . . . , A^ — 1 



WpAn) = ( ^ 



n 



yg-^'Xp_ — 



L yp - uJ^Xg 



Wp,q{n) = ifipfigT n 



ux, 
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where uj = e^'^*/^. Imposing ZTv-symmetry: Wp,q{n + A^) = Wp^g{n) (same for 
W), implies that Xq and yq are restricted to high genus Fermat-hke curves: 

M« + !) = < + <; A;xf = l-A:V-^; ky^ = l-k'^^^ (1) 

(the same also for Xp, yp, /ip). /c'^ = 1 — k"^ are temperature-like parameters, 
the low-temperature limit will be A;' — > . For given k, up to discrete choices 
of roots, eqs.(|l]) leave only two of the six functions independent, e.g. Xq and 

Xp. 

For N = 2 the chiral Potts model reduces to the Ising model and the Boltz- 
mann weights can be written in terms of Jacobi elliptic functions 



Xn 



vksn{q,k); i/q = vk cn{q,k)/dn{q,k), ^q = vk'/dn{q,k). 



W and W satisfy star-triangle-relations which involve the p- and g-variables 
separately, not as usual only through the difference p — q: 

N-l 

J]w,,,(/'-/)w^p,,(r'-/)Tl^p,,(/-rO=i?p,rW^p,,(/''-/0^p,r(/'-/''Ow^.,r(/'-/''')- 

Rpqr can be written in terms of the W, W and their iV-th roots [^ ^. From 
the Boltzmann weights, the transfer matrix is defined by 

L 

and the chiral Potts quantum chain hamiltonian arises as the derivative 

L N-l 



dTj 
q^p dq 



L- ^ ~ ^ = - E E ("' ^i + ^ "' ■'fi -^Tw') , (3) 



i=i 1=1 
where the operators Xj and Zj act in the spaces C^ at sites j and satisfy 



ZiX^ = X,ZiJ-^- 


Zf = Xf = 1; 


00 = e^-/^; 


(4) 


^i<p{2l-N)/N 


^i,t>{2l-N)/N 


A = l/k'. 


(5) 


smf 


smg 



and 

In right hand side of (^ we skipped an additive constant term. H commutes 
with the Z^r-charge operator Q = Y[j=i ^j- The eigenvalues of Q are u'^ with 
Q = 0, . . . ,N — 1. Since Tp ^ depends not only on p — q, it matters, at which p 
the limit g — *> p is taken. So, for given A, the hamiltonian (|^) depends on one 
chiral parameter 0. The coefficients (^ are conveniently written in terms of 
two chiral angles if and 0, which are related through 

COS(/9 = A COS0. (6) 



Using Lf as the independent variable, becomes complex for A < | cos(y9| and 
the hamiltonian (^ becomes non-hermit ian. The degrees of freedom remain 
unchanged, and the previously parity-nonconserving H now acquires a parity 
symmetry. Baxter mostly considers the 2-d model in this parameter range 
with (/? < |, since there the Boltzmann weights are positive real. In the single 
point A = 1 and </? = we recover the Fateev-Zamolodchikov parafermionic 
model pOi) see e.g. |H|. 



Performing a duality transformation on (^ we reach a model of the same type, 
which has and (/9 interchanged and A replaced by A~^ but also different b.c. 
Howes et al. [^ originally considered a self-dual model with ip = (p instead 



of (g). Probably this model is not integrable except at the line V' = f and the 
point v^ = 0, A = 1, where it coincides with (^. 

The analysis of the general 2-d model for A^ > 3 is difficult since for gen- 
eral k no simple hyperelliptic parametrization of the Wp^q, Wp^q is possible 



12|. Nevertheless, the free energy and the interface tension have been ob- 



tained circumventing an explicit uniformisation [0, ^. In [^ a gap-formula 



for the quantum chain has been obtained, from which the phase structure can 
be calculated: In the massive phases the ground state of the chain hamilto- 
nian is translational invariant (momentum P = 0). For the high-temperature 
(small A) massive phase the ground state is in the sector Q = 0. In the 
low-temperature phase varying A, ip or L, the ground state rotates among the 
P = 0, Q-sectors with many level-crossings. The incommensurate (IC) phase 
is reached when a P 7^ 0-state dips below the lowest P = state. The IC 
phase turns out to have a wedge shape, the tip of the wedge is at A = 1, (y? = 
(the parafermionic point) and spreading out with increasing (p. At </) = = | 
for A^ = 3 the IC-phase covers the interval 0.9013 < A < 1.1095. 

Much more detailed analytic information has been obtained for the special 
"superintegrable" case Xp = i/p, fip = I, or, in terms of the chain parameters: 
ip = (j) = ^. In this case H simplifies to 

j=i 1=1 
Splitting H^^'^ into two pieces by writing —H^^^ = Aq + XAi, one finds |^ that 
Aq and Al satisfy the Dolan-Grady |]19i conditions 



[Ao,[Ao,[Ao,Ai]]] = 16[Ao,Ai]; [A, [A, [Ai, Aq]]] = 16[A,Ao]. (8) 

Iff these conditions are met, starting from Ao, Ai one can generate Onsager's 
algebra A PT| : 



[Ai,Am\=4:Gi_m; [Gi, Am] = 2 Am+i - 2 Am-i; [Gi,Gm] = 0; l,meZ. 



A implies an infinite set of A-dependent commuting charges [Qm, Qm'] = : 

Qm = l{Am + A_m + HAm+i + A.rn+i)) among which Qo^H^'\ (9) 
As a consequence of A [|T8[ all eigenvalues E of H^^'' depend on A as 



niE 



E = a +bX + ^4mjVl + 2Acos0j + A2 (10) 

where a, b, rriE are integers, and one finds rrij = ±1/2. A gives no information 
about which values of a, b and 9j appear. Only few facts are known about the 



representation theory of A ||T8|, |T7 |. 



Starting from functional relations |jl|, |14|, |T^, and solving them via Bethe- 
ansatz-like equations Jl], |36|, the complete finite L spectrum for the Z3 su- 
perintegrable quantum chain has been given in [|16]. The general ZTv-case 
is discussed in p. All Z^v conserving boundary conditions (b.c.) X^+i = 
u^Xi] i? = 0, 1, . . . , A^— 1 can be treated. Eigenstates of H^'^^ are not known 



for general A. There are high- and low-temperature expansion results p4| 



3 Correlation functions 

To get more information about the phases of the Z^r-chiral Potts quantum 
chains we are interested in correlation functions. We assume periodic b.c. 
Xl+1 = Xi and will only consider expectation values in the lowest translational 
invariant state, which we denote by \v). This is the ground state in the non- 
IC regimes^. We concentrate on two-point functions of local operators Ej, 
where S^ stands for either Zj or Xj, or powers of these operators. We assume 
{v\v) = 1 and define as usual 

C^(r) = {v\E;^,E,\v) - {v\E^,\v) {v\E,\v) . (11) 

Since I f ) is translational invariant, C^{r) will be independent of j. We have 
Cxi-r) = Cxir)*; Cz{-r) = C|(r") = Cz(r). Thus, ahhough the chiral 
Potts model is not parity invariant, with periodic b.c. we can restrict ourselves 
to positive r with < r < L/2. According to the asymptotic behavior expected 
in the massive phases we will fit our numerical data at large r to the form 

C^ir) ^ as e-"/«= e^^'''^^^ + (1 - as)Sr,o- (12) 



^Albertini and McCoy Q have studied correlation functions in the superintegrable IC 
region via finite-size corrections to the energy levels and methods from conformal field theory. 



A finite A= means oscillatory correlations. In the high-temperature regime Cz 
is real and so it cannot show oscillations. Cz is real also in 8th-order low- 



temperature expansion |2^. So, we concentrate on Cx{t) and A will always 
mean Ax- Honecker [^ gave a rigorous argument that Cx must oscillate for 
(j) = 71 (this is a curve at A < 1 and </? > |) because of the shifted parity 
symmetry present there. Then, by continuity, we expect such oscillations in 
the adjacent part of the high-temperature phase too. For Z3 at = vr one 
gets |A| = 6, whereas from high-temperature expansion liniA^oA = Gn/^eip, 
so that on the superintegrable line for A — >■ we should have A ^ 12 and 
A = 00 at = 0. Observe that for small A a small change of ip produces a 
strong change of and so there we expect a fast 0-dependence of A. 
Information on A for A ^ 1+ can be obtained from a finite-size analysis of the 



pattern of ground-state level crossings. This way |21| concludes that in the 
thermodynamic limit A ~ (A — 1)~^/^ for A — »• 1_|_ and = |. 

Numerical calculations of Cx using Lanczos-type diagonalization of the 
sparce matrix H have been performed in |^, but since only chain lengths 
L < 14 or r < 7 can be managed, A (which, as we saw, should be 6... 12 or 
larger) could only be crudely determined. Here now the DMRG-method can 
do much better, and we will now turn to describe our new calculations. 



4 The DMRG-method 

For many years, say from 1980-1992, numerical investigations of quantum 
chain spectra mainly used Lanczos-type diagonalisation of finite chain hamil- 
tonians. In this method, the full space of states is considered, restricted only 
by possible conservation laws as e.g. momentum, Zjv-charge or parity. By an 
iterated application of the sparse hamiltonian matrix one achieves excellent 
convergence for typically about 10-20 lowest energy levels. For a Zs-model 
chains up to about L = 14 sites (ca.lO^ x 10^-matrices for periodic b.c.) can 
be treated without much effort. The extrapolation to the thermodynamic limit 
L —>■ 00 has to be based on these few L = 2, ... ,14 data. 
While in this method each chain length L is treated separately, the DMRG 



method invented by White [^] starts with small chains and enlarges the chain 
step by step. This method has become standard in condensed matter physics 
and we will be brief in summarizing the main idea at the specific example of 
a Za-chain with free b.c: 

We start with a small chain B^ of e.g. L = 5 sites. Then H is a mxm matrix, 
where m = 3^ = 243. Then we add one site and obtain a chain -B^+i of L -|- 1 



sites. We don't want to keep all 3m states of -Bl+i, but nevertheless want to 
get a good approximation for the ground state vector. 

In order to find out, which states of -Bl+i are most important for our goal 
(how to renormalise the system), we double the system by adding a system of 
the same size as shown in the Figure: 



Br 



Br 



Bl+1 
This 2(L + 1) sites system is called "superblock". We compute the 3^ m^ ~ 
5 ■ 10^-dimensional ground-state vector | $ ) of this superblock by a Lanczos 
algorithm: 

|$) = 5^<l>„-|2)^®|j)^ (13) 

and form the reduced density matrix taking the partial trace over the S-states: 

We diagonalise the reduced (729 x 729) density matrix pj^/ by some QL- 
algorithm and obtain the diagonal elements pi. Only those about m' ^ 80 
states with the largest pi are now kept to describe the renormalized matrix 
-Bi+i . This procedure is repeated until we reach chains of about L = 80 sites, 
still described by a m' x m' renormalised hamiltonian matrix. 
This method is most efficient for free b.c. If one wants to use periodic b.c, it 
is convenient to add a system repeated in the same orientation for not having 
two big adjacent blocks when joining the outer ends. We have not mentioned 
the decomposition into parity eigenstates, which is useful, even if parity is not 
conserved as it is in our chiral Potts model. For calculating matrix elements 
of operators, one has to keep track of the transformation of these operators 
into the new bases, which requires considerable effort. 

The DMRG method works very well even up to L = 80. ..200 sites keeping 
just about m' ~ 80 . . . 150 states. The reason for this is the exponentially de- 
creasing spectrum of the density matrix, which, as first pointed out by Nishino 



and Okunishi |30], can be understood from the close relation of the DMRG 
procedure to Baxter's variational method for calculating CTM eigenvalues [^]. 
We shall make use of this connection in the last section. 
Results for the correlation functions 

We have calculated Cx{x) in the high-temperature phase from chains up to 
L = 80 sites, so that it is reasonable to look up to r = 30 sites. The so-called 
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Figure 2: Real- and imaginary parts of the exponentially streched correlation 
function e^/^Cx(a;). Superintegrabel case </) = = |, high-temperature phase 
at A = 0.75. The fit gives A = 29.5(5) and i = 2.000(5). 



"infinite-size" DMRG approach |5^ is used, but with periodic b.c. (although 
this is much more effort in DMRG) in order to reduce boundary effects. Be- 
cause of space hmitations we give just one typical example for the superinte- 
grable case and A = 0.75. Of course, Cx decreases exponentially with r. So, in 
order to show the oscillation clearly, we multiply Cx by the factor e''/^ which 
can be done since our numerical precision is quite high. We show both the 
real and the imaginary part. 

We fit also 7 in ^^^ = 27 |1 — A| (in the superintegrable case the gap is linear 
in A). Our fits give 7 = l.OiO.l in lattice units. Going to smaller chiral angles 
(we give no Figures here, but these look similar to Fig. 2) we fitted the expo- 
nent z/ in ,^~^ ~ |1 — Al*^. z/ is found to change little with ip, not surprisingly 
since the mass-gap exponent ^ = § for the parafermionic case = is close 
to the value i^ = 1 at = |. However, A changes strongly: it must diverge for 
— > 0. In the integrable case (^ the behavior of A with decreasing is com- 
plicated because at A < 1 the range in A where the hamiltonian is hermitian, 
shrinks fast. However, for the Za-self-dual (and probably not integrable) case 
(j) = if within the errors we find A ~ 0~^. 



5 The Corner Transfer Matrix eigenvalues 

Baxter's corner transfer matrix (CTM) is a very efficient tool for the calcu- 
lation of single-spin expectation values (order parameters) of many integrable 
models. The CTM also plays also a crucial role in the Kyoto approach [^ to 
correlation functions of integrable models. 

A CTM is defined as the partition function of a corner of a planar lattice (e.g. 
quadrant in case of a square lattice, or a sextant in case of a triangular lat- 
tice [0). Fig. 3 shows a corner of a small square lattice. The empty circles 
indicate the spins forming the indices of the CTM, full circles are summed 
over, except for the lower diagonal spins, which are fixed or summed over, 
depending on whether we use fixed or free b.c. The CTM of the SE-corner of 




Figure 3: Corner transfer matrix A^s},{s'} 

a square lattice is usually called A, the other anticlockwise rotated ones B, C 
and D. The partition function is Z = Tr{ABCD). Baxter discovered that in 
the thermodynamic limit the spectrum of the CTM of the eight-vertex model 
takes a very simple form: let A^ be A diagonalised and normalized such that 
its maximum eigenvalue is unity. Then 

(Adjr^ri^) = exp {-Urx); rir e Z, (15) 

with X = 7iu/K and u is the rapidity parameter in the Boltzmann weights, K 
the elliptic half-period. For the derivation of (plSf ) the difference property of 
the Yang-Baxter equations is used and also some analyticity property which 
allows to make use of a periodicity property. From (pTS]) it is a short step to 
obtain the order parameter (magnetization). 

Of course, also for the Zn chiral Potts model one is interested in the order 
parameters. We like to know (Zq) for j = 1, . . . ,N — 1 where the spin with 
label is far from the boundaries. From series expansions and Z-invariance 



arguments there is the conjecture that the order parameters are independent 
of the chiral angle and their /c'-dependence is 

{zi) = {i~k'y^''-^'^'^^''"\ (16) 



Baxter |T0| , p!T|] has derived many important factorization and quasiperiodicity 
properties and functional relations of the chiral Potts CTM, but since the 
difference property is absent and analyticity seems to be restricted, as yet no 
exact proof of ( |l6|) was found. In the low-temperature limit Baxter did a 
truncated recursion calculation for the Za-case [0]. He obtains a 7th-order 
series approximation to the 17 largest eigenvalues Mi of M = ABCD. These 
eigenvalues depend only on a temperature-like parameter, arising from the 
nome of the hyperelliptic functions which parametrize the Boltzmann weights. 

"" M2''= Ms = r; M4 = M5 = 2r3 - Ir^ - -r^ - A^6 + IMI^r + . . . 

w T 2 27 ' 3 54 

M12 = (4 ± 2\/3)r^ + (10 ± 6\/3)r^ + . . . ; Mig = M17 = ^r^ + . . . . (17) 

15 2 

Here the largest eigenvalue is normalized to Mi = 1 and M2 is taken as the 
temperature scale r, which is related to k' of (|l]) or rather, u = k' /27, by 

r = m + 14m2 + 254m3 + 5150m^ + 111123u^ + 2495269m^ + 57608712u^+. . . (18) 

In contrast to other integrable models with difference property, here the spec- 
trum of log M is not equidistant. Baxter uses fixed spin zero b.c. on the CTM, 
so the ground state is non-degenerate. The degenerate levels are those with 
center spin 1 and 2, the non-degenerate ones have center spin zero. Baxter's 
derivation of (|l3) is best justified in the region of real positive Boltzmann 
weights and it is not clear which details can be trusted when extending the 
formulae to the region of complex Boltzmann weights. Using (|T7|), ([18|) one 
finds correctly (Zg) = k' ' up to the order r^ considered, in agreement with 
([T6| ) as first conjectured from the superintegrable chain. 

Looking into the apparent convergence of the expansions ([T7|), (p!8D, one is 
inclined to trust the lowest eigenvalues M2 . . . M5 up to about k' ^ 0.7 within 
a few percent: For k' = 0.7 one finds M2 = t ^ 0.025. However, approaching 
A;' = 1 we clearly see that convergence of the series gets lost quite suddenly, 
although in 7th-order at k' = 0.9 we still remain at t ^ 0.06. 



6 CTM eigenvalues from DMRG density matrices 



Fig. 4 illustrates that the matrix M = ABCD is just the partition function 
of a lattice with a seam opened between A and D. There is a direct relation 




Figure 4: The product of four corner transfer matrixes {ABCD)i^„^,{rj'}- The 
seam in the lower vertical indicates that no trace is taken. 

between M and the reduced density matrix pj j/ of (p^): Baxter 0, chapt.13, 
pointed out that in the thermodynamic limit the product of two CTM, say 
AB in Fig.4, is the maximal eigenvector | \E' ) of the column-to-column transfer 
matrix (TM) (here moving left to right). Now, for integrable models the TM 
of the 2-dim. system is in one-to-one correspondence with the quantum chain 
hamiltonian H as e.g. in (j^) (vice-versa one uses the Trotter construction) 
and the maximal eigenvector of the TM is the ground state eigenvector of if, 
i.e. I \E' ) ~ I $ ) of eq.(|13|). So M = ABCD, which is summed over the closed 
seam betwen B and C, is proportional to pj j/ of eq.(0). All this is correct for 
large systems, i.e. when the boundaries play no essential role. 

We now compare numerically the pi^i' (which appeared as a technical tool 
for selecting the most relevant states in our DMRG calculation) to Baxter's 
CTM eigenvalues (p!7[). Baxter uses fixed b.c. at low k' . It is not easy to 
implement fixed b.c. in DMRG, so we use free b.c. for the dual model P] at 
high temperatures A < 1. We have studied several values of L and y?, A. In 
Fig. 5 we show that there is excellent agreement already when we use a small 
L = 10 system, which can be understood considering that for A < 0.75 the 
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Figure 5: The largest 15 eigenvalues for several A. Discrete points: Baxter's 
CTM-calculation (p^ together with our numerical density matrix eigenvalues 
for free b.c. and L = 10: solid lines: at the boundary of the hermitian region 
= 0, ip = arccos A; dashed lines: on the superintegrable line. 



correlation length ^ is smaller than three lattice sites. E.g. for A = 0.5 eqs.(]T7|) 
give M4 = 2.44884-10"^ whereas the diagonalisation of the density matrix gives 
P4 = 2.4418 • 10"^ at = f and p4 = 2.4382 • 10"^ at = 0. Baxter's formulae 
are not valid at A = 1, so any agreement there is not significant. For the 
Ising and XXZ-models, the checks analogous to the ones given here were made 
earlier by Peschel et al. ^^. In their cases the CTM-eigenvalues are known 
exactly and, apart from degeneracies, are strictly logarithmically equidistant. 

G.v.G. gratefully acknowledges the support by INTAS-97-1312. 
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